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EXECUTIVE  SUMMARY 


OBJECTIVE 

The  objective  of  this  study  was  to  determine  the  grazing  optical  reflectivity  of  an  ocean  covered  by 
capillary  waves. 

APPROACH 

Starting  with  an  integral  equation  for  ocean  reflectivity,  we  introduced  approximations  applicable 
in  the  grazing  regime  and  arrived  at  an  analytic  formula  for  the  grazing  ocean  reflectivity  of  circular 
sources  such  as  the  Sun. 

RESULT 

The  analytic  approximation  depends  on  the  variance  of  the  capillary  wave  slopes,  the  angular 
radius  of  the  source,  and  the  grazing  angles  of  the  source  and  receiver. 


% 


ili 


Preceding  Pagcfelank 


CONTENTS 


EXECUTIVE  SUMMARY . '» 

1.  INTRODUCTION . 1 

2.  OCEAN  REFLECTIVITY . 3 

3.  GRAZING  OCEAN  REFLECTIVITY . 5 

4.  SMOOTH  OCEAN  REFLECTIVITY . 7 

5.  ROUGH  OCEAN  REFLECTIVITY . 13 

6.  FINAL  RESULT . 17 

7.  REFERENCES . 21 

APPENDIX  A:  GRAZING  AREA  OF  A  ROUGH  SURFACE . . A-1 

Figures 

1 .  Grazing  geometry  for  capillary  wave  reflectivity.  A  specular  reflection  occurs  at  the  flat 
capillary  wave  facet  that  is  tilted  with  respect  to  the  zenith  direction  by  the  angle,  0n. 

The  source,  assumed  to  be  circular,  has  an  angular  radius  of  e.  The  facet  is  shown 

with  a  positive  slope  in  the  Y  direction . 5 

2.  Plass-Cox-Munk  interaction  probability  density,  qr ,  for  a  wind  speed  of  10  m  s'1  and  a 
receiver  elevated  by  4.4  mrad  (1/4°).  The  receiver  looks  directly  toward  the  Sun.  The 
figure  is  drawn  for  a  smooth  surface  reflectivity  of  unity.  Each  dark  area  represents  the 
top  of  a  sun  glint  column  for  solar  elevations  increasing  from  0  to  210  mrad  (0°  to  12°) 

in  equal  steps  of  35  mrad  (2°) . 6 

3.  Smooth  surface  reflectivity  of  seawater  as  a  function  of  average  elevation  for  an  optical 
wavelength  of  3.2  pm.  Circles  are  the  result  of  an  exact  calculation  and  lines  show  an 
exponential  fit  given  by  equation  (7).  The  data  labeled  “S”  are  for  radiance  polarized 
perpendicular  to  the  plane  of  incidence,  the  data  labeled  “P”  are  for  radiance  polarized 

parallel  to  the  plane  of  incidence,  and  the  unlabeled  data  are  for  unpolarized  radiance . 8 

4.  Smooth  surface  reflectivity  of  seawater  as  a  function  of  optical  wavelength  for  an  average 

elevation  of  100  mrad.  The  curves  have  the  same  meaning  as  those  given  in  the  previous  figure. 
The  vertical  dashed  line  indicates  the  wavelength  used  in  figure  3 . 9 

5.  A  sequence  of  tolerance  ellipses  for  the  rising  (or  setting)  Sun.  The  receiver  is  looking 
directly  toward  the  Sun  from  an  elevation  of  1  mrad.  The  parameter  next  to  each  ellipse 
is  the  elevation  of  the  solar  center  in  terms  of  the  solar  radius,  e.  The  dashed  line  shows 

the  approximation  for  the  semi-major  axis  of  the  ellipse  given  by  equation  (11) . 1 5 

6.  Rough  ocean  reflectivity  as  a  function  of  solar  elevation.  The  receiver  is  elevated  by 
1  mrad  and  the  wind  speed  is  10  m  s'1.  The  solid  line  labeled  “Exact”  is  a  numerical 
calculation  based  on  equations  (1)  and  (2).  The  dashed  line  labeled  “Rough”  is 
equation  (14).  The  dashed  line  labeled  “Flat-Top"  is  equation  (23)  of  Zeisse  (1995). 

The  dashed  line  labeled  “Cox-Munk/50”  is  equation  (9)  of  Cox  and  Munk  (1954)  after 
division  by  a  factor  of  50 . 18 

7.  The  same  as  figure  6  for  a  wind  speed  of  1/4  m  s'1 . 19 

v 


Tables 

Fresnel  reflectivity  of  seawater  at  an  average  grazing  angle  of  100  mrad.  The  constants, 
bs  and  bp,  are  defined  in  equation  (7) . 


10 


vi 


1.  INTRODUCTION 


The  reflection  of  distant  objects  in  the  sky  by  the  rough  sea  surface  produces  interesting  natural 
effects.  Several  thoughtful  observers  have  commented  on  these  effects.  Minnaert  (1954)  says,  “The 
first  25°  or  35°  of  the  sky  above  the  horizon  are,  in  reality,  hardly  visible  in  the  reflection.  [This 
phenomenon  is]  easily  explained;  at  a  great  distance  we  can  only  see  the  sides  of  the  waves  turned 
towards  us.  This  makes  it  seem  as  if  we  saw  all  the  objects  in  the  sky  reflected  in  a  slanting  mirror.” 
Cox  and  Munk  (1955)  pointed  out  that  “If  the  sea  were  absolutely  flat,  then  the  radiance  of  the  sea 
surface  just  beneath  the  horizon  would  equal  the  radiance  of  the  sky  just  above  it,  and  there  would  be 
no  visible  horizon.”  In  other  words,  the  fact  that  we  see  the  ocean  horizon  is  proof  that  the  sea  is 
rough.  Shaw  (1999)  recently  re-emphasized  that  the  width  of  a  sun  glint  pattern  is  a  sensitive  measure 
of  the  wind  speed  directly  above  that  pattern.  All  these  comments  refer  to  a  grazing  geometry  when 
the  source  and  receiver  are  at  relatively  low  altitudes,  say  several  tens  of  meters,  but  are  separated  by 
a  range  of  several  tens  of  kilometers.  This  report  presents  an  approximate  formula  for  the  grazing 
reflection  of  a  circular  source  by  the  rough  sea.  The  formula  can  be  applied,  for  example,  to  the 
reflectivity  along  the  center  of  a  sun  glint  pattern  or  to  the  reflectivity  of  a  missile  plume  (considered 
as  a  circle)  as  seen  from  the  deck  of  a  ship. 

In  the  optical  region,  the  reflectivity  of  an  ocean  without  gravity  waves  can  be  calculated  using 
geometrical  optics.  The  reason  is  that  the  wavelength  of  the  radiation  being  reflected  is  much  less 
than  the  radius  of  curvature  of  a  typical  capillary  wave  facet.  The  ocean  can  therefore  be  modeled  as  a 
collection  of  flat,  mirror-like  facets  whose  tilts  fluctuate  at  random  under  the  wind’ s  influence.  If  the 
slope  of  a  facet  is  correct,  a  ray  from  a  source  such  as  the  Sun  will  be  reflected  into  an  optical  receiver 
such  as  the  eye.  The  problem  then  reduces  to  finding  the  Fresnel  reflectivity  of  a  small  seawater 
mirror  and  counting  the  relative  number  of  mirrors  which,  on  average,  can  reflect  a  ray  from  the 
source  into  the  receiver. 

Of  course,  Cox  and  Munk  (1954,  1955,  and  1956)  have  already  addressed  this  problem  in  a 
famous  series  of  papers.  They  measured  the  statistical  occurrence  of  various  capillary  wave  slopes  by 
taking  photographs  of  sun  glitter  from  an  airplane.  However,  their  equation  for  reflected  radiance 
applies  to  the  case  where  the  ocean  is  observed  from  above,  not  to  the  grazing  case.  When  calculating 
reflectivity,  Cox  and  Munk  assumed  that  a  patch  of  ocean  of  area  A  appears  to  have  an  area  of  A  sin'F 
when  viewed  from  the  side  at  a  grazing  angle,  VF.  In  other  words,  for  that  part  of  their  calculation, 
they  assumed  that  the  patch  was  flat.  Because  radiance  involves  the  power  received  per  unit  solid 
angle  per  unit  projected  area,  and  because  the  projected  area  of  an  ideal  flat  surface  approaches  zero 
at  a  grazing  angle  of  observation,  the  radiance  obtained  by  Cox  and  Munk  approaches  infinity  at 
grazing.  Even  so,  their  analysis  is  quite  good  for  those  elevations  greater  than  about  10°.  For  observa¬ 
tions  at  elevations  less  than  10°,  especially  for  those  approaching  0°,  their  approach  must  be  modified 
by  using  the  fact  that  the  sea  is  a  rough  surface  that  can  be  seen  even  at  very  small  grazing  angles. 

It  is  interesting  to  compare  the  results  already  obtained  by  the  radar  community  for  the  grazing 
reflectivity  of  the  rough  ocean  surface.  Radar  wavelengths  fall  in  the  range  of  1  mm  to  1  cm,  and  a 
typical  radar  footprint  (Kerr,  Fishback,  and  Goldstein,  1951)  (given  by  the  size  of  the  first  Fresnel 
zone)  is  about  1  km  along  the  propagation  direction  and  10  m  perpendicular  to  it.  The  basic  radar 
reflection  model  (Ament,  1953)  assumes  that  the  ocean  is  covered  by  gravity  waves  alone.  With  this 
assumption,  geometrical  optics  can  be  used  for  radar  because  the  radar  wavelength  is  much  less  than 
the  wavelength  of  an  ocean  gravity  wave.  The  sinusoidal  gravity  wave  is  approximated  by  a  series  of 
horizontal  sections  and  a  specular  reflection  is  assumed  to  occur  at  each  section.  Each  horizontal 
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section  with  altitude,  y,  introduces  a  phase  shift,  Ak  y  sin  /X ,  compared  to  the  reflection,  F0,  from 
the  smooth  surface  at  y  =  0  (mean  sea  level).  After  superimposing  gravity  waves  at  various  wave¬ 
lengths,  the  sections  are  assumed  to  have  an  altitude  that  is  randomly  distributed  according  to  a 
Gaussian  law  about  mean  sea  level.  The  total  rough  surface  reflectivity,  T,  is  found  by  integrating 
over  all  possible  altitudes  with  the  result,  for  radar,  that 

p 

—  =  p  =  exp(-g2/2), 

(1) 

g  =  An  <rh  sin  *¥/ X, 

where  Oh  is  the  standard  deviation  of  the  gravity  wave  altitude  distribution.  This  result  underestimated 
measurements  of  very  rough  surfaces  and  was  improved  by  assuming  that  the  gravity  waves  not  only 
had  a  Gaussian  amplitude  distribution,  but  also  a  uniform  phase  distribution  (Miller,  Brown,  and 
Vegh,  1984). 

These  considerations  for  radar  apply  only  to  so-called  coherent  (specular)  reflection;  there  is  also  a 
diffusely  reflected  (scattered)  component.  A  recent  review  of  radar  rough  surface  reflectivity  (Barrick, 
1998)  notes  that  propagation  and  scattering  are  “inextricably  connected”  at  grazing  and  that  their 
improper  separation  has  led  to  contradictions  in  associating  theory  with  experiment. 

From  my  point  of  view,  the  primary  distinction  between  radar  reflectivity  and  optical  reflectivity  is 
this:  For  radar  reflectivity  one  assumes  that  local  regions  of  the  reflecting  surface  are  horizontal 
whereas  for  optical  reflectivity  one  assumes  that  local  regions  are  tilted.  When  coupled  with  the  fact 
that  the  Fresnel  reflectivity  of  any  smooth  surface  approaches  100%  at  grazing,  this  results  in 
qualitatively  different  behavior  for  grazing  reflectivity  of  rough  radar  and  rough  optical  surfaces.  In 
the  most  coarse  approximation,  the  radar  sea  is  nominally  flat  and  the  reflectivity  at  exact  grazing 
(4*  =  0)  is  close  to  100%  whereas  the  optical  sea  is  mostly  tilted  and  the  reflectivity  at  exact  grazing  is 
close  to  0%. 
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2.  OCEAN  REFLECTIVITY 


Projection  of  a  rough  surface  toward  the  observer  and  the  neglect  of  gravity  waves  are  important 
aspects  of  the  optical  description  of  ocean  reflectivity.  The  second  aspect  is  a  severe  restriction. 
Gravity  waves  significantly  influence  the  tilt  of  all  capillary  wave  facets  and,  hence,  significantly 
influence  the  reflectivity.  Furthermore,  I  do  not  consider  shadowing  or  multiple  reflections,  two 
effects  that  become  increasingly  important  as  exact  grazing  is  approached.  Therefore,  I  will  be 
content  with  a  rough  form  for  grazing  reflectivity,  perhaps  good  to  within  only  a  factor  of  two  for 
larger  sources  such  as  the  Sun  at  lower  wind  speeds  such  as  lA  m  s'1. 

Let  T  be  the  optical  reflectivity  of  the  rough  ocean  surface.  This  means  that  when  a  radiance,  Na,  is 
incident  on  the  ocean  surface,  the  reflected  radiance,  Nr,  will  be  given  by  Y  Na.  The  following  pair  of 
equations  (Zeisse,  1995)  gives  a  rough  surface  reflectivity  that  is  the  mathematically  exact  conse¬ 
quence  of  the  physical  approximations  mentioned  in  the  previous  paragraph: 

r=^-=  JJ  r.  ,r  (i) 

V  o  ellipse 

U  r-const 

(cos  (O  cos  0n  )p 

<J r  —  - \  *  \^/ 

JJ  (cos m  cosen)pdCxd^y 

0)<JT  2 

Ur=const 


In  these  equations,  co  is  the  angle  of  incidence  for  a  specular  reflection  between  source  and  receiver  at 
a  capillary  wave  facet,  T0  (co,  A)  is  the  smooth  surface  reflectivity  of  seawater  at  co  and  optical 
wavelength,  A,  0n  is  the  zenith  angle  of  the  facet  normal  (also  known  as  the  facet  tilt),  and  C,x  and  £y 
are  the  facet  slopes  in  the  upwind  and  crosswind  directions,  respectively.  The  function,  p( £y>  W), 
is  the  well-known  Cox-Munk  wave  slope  probability  density. 
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-exp 


_L 

-2 


<?l  =  0.000  +  3.16x10-3W 
crc2  =  0.003  +  1. 92  x10-31¥ 


(3) 


which  has  the  shape  of  a  two-dimensional  bell  in  slope  space.  In  equation  (3),  <7„2  and  <rc2  are  the  facet 
slope  variances  in  the  upwind  and  crosswind  directions,  respectively,  and  W  is  the  wind  speed  in  m  s 
\  The  quantity  W)  d£x  dC,y  gives  the  probability  that  a  wave  facet  will  have  a  slope  within  + 

d£t/2  of  £  and  +  d^/2  of  £y  when  the  wind  speed  is  W.  This  quantity  can  also  be  interpreted  as  the 
fractional  area  within  the  footprint  occupied  by  all  facets  whose  slope  falls  within  the  specified  range. 


The  function  ^embodies  the  assumption,  first  put  forward  by  Plass,  Kattawar,  and  Guinn  (1975), 
that  the  importance  of  a  particular  facet  within  a  wind-ruffled  ocean  footprint  is  given  by  the  area  of 
that  facet  projected  toward  the  receiver.  For  a  single  facet  whose  slope  is  {^x,  £  v ),  this  projection  is 

accomplished  by  the  factor  (cos®/cos0„  )  in  equation  (2).  Therefore,  the  denominator  of  qr 

represents  the  area  of  all  facets  projected  toward  the  receiver,  which  is  the  same  as  the  projected  area 
of  the  entire  footprint  considered  as  a  rough  surface.  When  performing  the  integral  in  the 
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denominator  of  qr,  the  notations  under  the  integral  signs  serve  as  a  reminder  that  (1)  the  angle  of 
incidence  must  remain  less  than  90°,  and  (2)  receiver  coordinates  are  held  fixed  during  the 
integration.  The  first  restriction  means  that  only  the  front  of  each  facet  is  used  to  evaluate  the  integral 
because,  as  Minnaert  (1954)  says,  “at  a  great  distance  we  can  only  see  the  sides  of  the  waves  turned 
towards  us.” 

Finally,  the  word  “ellipse”  under  the  integral  in  equation  (1)  means  that  the  integral  is  confined  to 
slopes  within  the  tolerance  ellipse1.  The  area  of  the  tolerance  ellipse  is 

E=  JJ 

ellipse 

where  £  is  the  angular  radius  of  the  source  subtended  at  the  ocean  footprint.  The  approximation  made 
in  equation  (4)  is  that  the  quantities  inside  the  bracket,  namely,  the  angle  of  incidence  and  the  facet 
tilt,  remain  constant  as  the  ray  from  the  source  varies  across  the  disk. 

Equations  (1)  and  (2)  apply  at  all  receiver  elevations,  but  they  are  obscure.  Fortunately,  at  receiver 
elevations  greater  than  about  10°  they  reduce  to  the  simple  form  used  by  Cox  and  Munk  (1954, 1955, 
and  1956).  In  the  grazing  regime  equations  (1)  and  (2)  must  still  be  used  and  there  they  must  be 
numerically  evaluated.  In  this  report,  I  derive  a  more  easily  evaluated  approximate  form. 


sec  (O  sec3  &„ 


>7t£ 


(4) 


1  When  there  is  a  specular  reflection  at  the  facet  from  the  sky  to  a  receiver  whose  position  is  held  fixed,  the  law  of 
reflection  can  be  used  to  associate  a  sky  location  with  an  ocean  slope.  By  considering  the  Jacobian  of  the 
transformation  from  ocean  slopes  to  a  coordinate  system  centered  on  a  circular  source  in  the  sky,  such  as  the  Sun,  it 
can  be  shown  that  a  disk  in  the  sky  is  associated  with  an  elliptical  region  of  slope  space.  Cox  and  Munk  called  this 
region  the  “tolerance  ellipse.” 
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3.  GRAZING  OCEAN  REFLECTIVITY 


Consider  a  receiver  (such  as  a  camera  in  the  infrared  or  an  eye  in  the  visible)  observing  a  distant 
circular  source  (such  as  the  Sun  or  the  mirror  of  a  telescope  transmitter).  Imagine  that  the  camera  is 
mounted  on  a  ship,  or  the  observer  is  on  the  shore,  while  the  source  is  near  the  horizon  in  the  distant 
sky.  Figure  1  shows  the  geometry  for  this  situation,  where  4^  and  4*r  are  the  source  and  receiver 
elevations,  respectively.  I  chose  a  coordinate  system2  that  has  an  origin  at  the  point  of  reflection,  with 
the  Z-axis  pointing  to  the  zenith,  the  Y-axis  pointing  to  the  center  of  the  source,  and  the  X-axis  such 
that  a  right-handed  system  is  formed.  Let  U  be  a  unit  vector  with  coordinates  (a,  b,  c)  in  this  system. 


2  e 
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Figure  1 .  Grazing  geometry  for  capillary  wave  reflectivity.  A  specular  reflection  occurs  at  the  flat 
capillary  wave  facet  that  is  tilted  with  respect  to  the  zenith  direction  by  the  angle,  0n.  The  source, 
assumed  to  be  circular,  has  an  angular  radius  of  e.  The  facet  is  shown  with  a  positive  slope  in  the 
Y  direction. 


The  law  of  reflection  provides  the  following  relationships  among  the  angles  in  figure  1 : 
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(5) 


where  4*  is  the  average  grazing  angle  of  the  source  and  receiver.  Typical  values  for  grazing  angles 
can  be  obtained  from  the  following  considerations.  The  height  of  the  receiver  may  vary  from  several 
tens  of  meters  (for  a  camera  mounted  on  a  ship)  to  several  meters  (for  a  human).  In  either  case,  the 
range  from  the  receiver  to  the  ocean  footprint  will  vary  from  several  km  to  several  tens  of  km.  The 
elevation  of  the  center  of  the  source  may  vary  from  a  minimum  when  the  limb  of  the  source  just 
touches  the  geometrical  horizon  to  an  arbitrary  maximum  of  about  100  mrad  (5.73°).  Hence,  the 
grazing  angles  will  fall  in  the  following  ranges  for  the  case  when  the  source  is  the  Sun: 

1  mrad  <  4L  <10  mrad 

(< 

4  mrad  <  4^  <  100  mrad 


2  In  previous  equations,  here,  and  in  Zeisse  (1995)  and  Cox  and  Munk  (1954, 1956),  the  X-axis  pointed  in  the 
upwind  direction. 
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From  this  point  on,  I  assume  that  the  receiver  is  looking  directly  toward  the  source,  as  figure  1  shows, 
and  I  consider  grazing  angles  in  the  ranges  given  by  equation  (6).  I  also  assume  that  the  source  is  the 
Sun,  although  the  derivation  applies  to  a  circular  source  of  any  radius  up  to  and  including  the  Sun’s 
radius  (4.641  mrad). 

Figure  2,  reprinted  from  Zeisse  (1995),  shows  how  qr  depends  on  slope  when  the  wind  speed  is 
10  ms"1  and  the  ray  propagating  towards  the  receiver  leaves  the  footprint  at  the  grazing  angle  of  14°. 
The  dark  regions  shown  on  the  surface  of  qr  are  the  tops  of  a  parade  of  sun  glint  columns  for  seven 
different  positions  of  the  Sun.  A  sun  glint  column  is  that  portion  of  slope  space  over  the  tolerance 
ellipse  and  under  the  function,  T0  qr.  The  entire  volume  under  ro  qr  is  unity  according  to  equation  (2), 
and  the  volume  of  a  glint  column  is,  according  to  equation  (1),  the  rough  surface  reflectivity  I  seek. 


Figure  2.  Plass-Cox-Munk  interaction  probability  density,  qr,  for  a  wind  speed  of  10  m  s"1  and  a 
receiver  elevated  by  4.4  mrad  (1/4°).  The  receiver  looks  directly  toward  the  Sun.  The  figure  is  drawn 
for  a  smooth  surface  reflectivity  of  unity.  Each  dark  area  represents  the  top  of  a  sun  glint  column  for 
solar  elevations  increasing  from  0  to  210  mrad  (0°  to  12°)  in  equal  steps  of  35  mrad  (2°). 
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4.  SMOOTH  OCEAN  REFLECTIVITY 


Figure  3  shows  the  smooth  surface  reflectivity  of  seawater.  The  unpolarized  values  are  the  average 
of  the  two  polarized  values,  “S”  for  radiance  polarized  normal  to  the  plane  of  incidence,  and  “P”  for 
radiance  polarized  parallel  to  the  plane  of  incidence.  The  circles  are  the  result  of  an  exact  calculation 
(Stratton,  1941)  for  the  Fresnel  reflectivity  of  a  medium  with  the  complex  optical  index  (Hale  and 
Querry,  19733;  Querry  et  al.,  19774)  of  seawater.  The  lines  for  the  polarized  data  in  figure  3  are  an 
exponential  fit  that  is  forced  to  go  between  1  at  an  elevation  of  zero  and  the  exact  value  at  an 
elevation  of  100  mrad.  The  line  for  the  unpolarized  data  is  the  average  of  the  other  two  lines.  That  is, 

Ros  0?,  A)=exp[-^(A>P] 

^)=exp[-&p(A)4/] 

A)=  foOF,  AktfJ'F,  A)}/2  (7) 

b,=-  lOxlnrJO.U) 

bp  = -ioxinrop(o.i,  A) 

where  ro  ('P,  A)  is  the  calculated  value  and  R0  (41,  A)  is  the  fitted  value.  (I  have  used  equation  (5)  to 
replace  the  angle  of  incidence  with  the  grazing  angle  in  these  arguments.)  For  wavelengths  between 
0.4  and  12  pm,  and  elevations  between  0  and  100  mrad,  the  maximum  absolute  error  between  the 
calculation  and  the  fit  is  0.0016.  This  maximum  error  occurs  at  the  wavelength  of  3.2  pm  used  for 
figure  3.  Figure  4  shows  how  the  Fresnel  reflectivity  depends  on  wavelength  and  table  1  lists  the 
value  of  the  constants,  bs  and  bp,  as  a  function  of  wavelength. 


3  These  data,  for  pure  water,  were  used  from  0.2  to  2  pm. 

4  These  data,  for  seawater,  were  used  from  2.2  to  12  pm. 
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AVERAGE  ELEVATION  *F  (mrad) 

Figure  3.  Smooth  surface  reflectivity  of  seawater  as  a  function  of  average  elevation  for  an  optical 
wavelength  of  3.2  pm.  Circles  are  the  result  of  an  exact  calculation  and  lines  show  an  exponential  fit 
given  by  equation  (7).  The  data  labeled  “S”  are  for  radiance  polarized  perpendicular  to  the  plane  of 
incidence,  the  data  labeled  “P”  are  for  radiance  polarized  parallel  to  the  plane  of  incidence,  and  the 
unlabeled  data  are  for  unpolarized  radiance. 
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FRESNEL  REFLECTIVITY  r 


WAVELENGTH  X  (tmi) 

Figure  4.  Smooth  surface  reflectivity  of  seawater  as  a  function  of  optical  wavelength  for  an  average 
elevation  of  100  mrad.  The  curves  have  the  same  meaning  as  those  given  in  the  previous  figure.  The 
vertical  dashed  line  indicates  the  wavelength  used  in  figure  3. 
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Table  1 .  Fresnel  reflectivity  of  seawater  at  an  average  grazing 
angle  of  100  mrad.  The  constants,  bs  and  bp,  are  defined  in 
equation  (7). 


A,  (|im) 

bs 

bp 

0.4 

4.48 

8.10 

0 . 6 

4.53 

8.11 

0.8 

4.55 

8.11 

1.0 

4.55 

8.11 

1.2 

4.59 

8.12 

1.4 

4.62 

8.13 

1.6 

4.65 

8.14 

1.8 

4.69 

8.15 

2.0 

4.74 

8.16 

2.2 

4.77 

8.17 

2.4 

4.92 

8.22 

2.6 

5.30 

8.36 

2.8 

5.86 

8.31 

3.0 

3.54 

7.56 

3.2 

3.61 

8.05 

3.4 

3.95 

8.05 

3.6 

4.14 

8.06 

3.8 

4.27 

8.07 

4.0 

4.36 

8.08 

4.2 

4.41 

8.09 

4.4 

4.47 

8.10 

4.6 

4.47 

8.10 

4.8 

4.50 

8.10 

5.0 

4.53 

8.11 

5.2 

4.59 

8.12 

5.4 

4.69 

8.15 

5.6 

4.83 

8.19 

5.8 

5.08 

8.26 

6.0 

4.78 

8.02 

6.2 

4.16 

7.99 

6.4 

4.34 

8.06 

6 . 6 

4.45 

8.08 

6.8 

4.53 

8.10 

7.0 

4.57 

8.11 

7.2 

4.61 

8.12 

7.4 

4.66 

8.13 

7.6 

4.71 

8.14 

7.8 

4.76 

8.16 

8.0 

4.80 

8.17 

8.2 

4.85 

8.18 

8.4 

4.90 

8.19 

8.6 

4.96 

8.21 

8.8 

5.02 

8.23 

9.0 

5.05 

8.24 

9.2 

5.09 

8.25 

9.4 

5.15 

8.27 

9.6 

5.25 

8.31 

9.8 

5.34 

8.35 
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Table  1 .  Fresnel  reflectivity  of  seawater  at  an  average  grazing 
angle  of  100  mrad.  The  constants,  bs  and  bp,  are  defined  in 
equation  (7).  (continued) 


X  ( |Ltm) 

bs 

bp 

10.0 

5.47 

8.40 

10.2 

5.59 

8.44 

10.4 

5.74 

8.50 

10.6 

5.88 

8.55 

10.8 

5.97 

8.55 

11.0 

6.00 

8.51 

11.2 

5.96 

8.40 

11.4 

5.80 

8.20 

11.6 

5.55 

7.93 

11.8 

5.26 

7.64 

12.0 

4.92 

7.32 
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5.  ROUGH  OCEAN  REFLECTIVITY 


In  this  section  (and  the  appendix),  I  introduce  approximations  into  equations  (1)  and  (2)  that  reduce 
them  to  an  easily  evaluated  algebraic  form.  The  final  approximate  form,  equation  (14),  depends  only 
on  the  grazing  angles,  the  source  radius,  and  the  wind  speed. 

The  first  simplification  eliminates  the  distinction  between  upwind  and  crosswind  slope  variances. 

I  replace  each  of  these  with  their  arithmetic  mean  value: 


a 


2 


2 

0.0015  +  2.54X10-3  W. 


(8) 


One  standard  deviation,  containing  90%  of  all  wave  slopes,  corresponds  to  a  tilt  of  6n  =  tan  1 G  ~ 

64  mrad  (3  2/3  °),  at  a  wind  speed  of  1  m  s'1  or  163  mrad  (9  1/3  °)  at  a  wind  speed  of  10  m  s'1. 

After  making  this  change  and  inserting  equations  (2)  and  (3)  into  equation  (1),  the  reflectivity  is 

r  =  — — —  JJ  T0  cos  (£>  sec0„  exp{-  (i u 2  +  v2  )/2 }du  dv ,  (9) 

ellipse 

Ur=const 


where  the  variables  of  integration  are 


V  =  Cy/0- 


(10) 


I  have  also  replaced  the  denominator  of  equation  (2)  by  its  geometrical  representation,  BJA,  where  A 
is  the  footprint  area  in  the  plane  of  mean  sea  level  and  B,  i s  the  area  (projected  towards  the  receiver) 
of  the  rough  sea  surface  above  that  footprint.  During  the  integration  of  equation  (9),  the  receiver  ray 
stays  fixed  while  the  source  ray  varies  systematically  across  the  Sun  with  the  facet  tilt  adjusting  itself 
to  maintain  a  specular  reflection. 

To  understand  the  behavior  of  equation  (9),  we  must  understand  how  the  tolerance  ellipse  appears 
in  slope  space  during  grazing.  Figure  5  shows  how  the  shape  of  the  tolerance  ellipse  depends  on  solar 
elevation  when  the  receiver  is  elevated  by  1  mrad.  An  estimate  of  a,  the  semi-major  axis  of  the 
tolerance  ellipse,  may  be  obtained  from  an  estimate  of  the  semi-minor  axis  b  and  approximate  area  of 
the  ellipse  given  by  equation  (4): 


-  x  x  j 

secosec  0„£_  e 
2 

When  the  lower  limb  of  the  Sun  just  touches  the  geometrical  horizon,  and  for  the  receiver  at  exact 
grazing  0Fr  =  0),  Y  =  TV  2  =  e  /  2  and  equation  (11)  shows  that  a  ~  1,  corresponding  to  a  tilt  of  45°. 
At  grazing,  the  tolerance  ellipse  has  become  very  large  in  the  X-direction,  so  large  that  the  trigono¬ 
metric  terms  in  equation  (9)  vary  considerably  from  their  value  at  the  center  of  the  ellipse.  Never¬ 
theless,  I  may  remove  these  terms,  and  T0  as  well,  from  the  integral  in  equation  (9)  without 
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appreciable  error  because  these  terms  vary  slowly  across  the  ellipse  compared  to  the  rapid  fall-off  of 
the  exponential  term  away  from  the  center.  After  removing  them,  replacing  them  by  their  central 
value,  and  then  applying  the  small  angle  approximation  to  the  central  value,  equation  (9)  becomes 

r  »  —  JJ  exp{-  (m 2  +  v2  )/ l\du  dv .  (12) 

ellipse 

Ur-const 


Next,  I  replace  the  ellipse  by  the  rectangle  that  just  contains  it  and  introduce  a  factor  of  7t/4  to  take 
into  account  the  different  areas  of  these  two  shapes: 


r  = 


Br  2 n  4 


+b  a 


Jrf«exp(-M2/2)  Jdvexp(-v2/2) 

cr 


-b  cr 


A  r-*2«  erfU-  W 


B ,  8 


4lc  J J  {420  J 


«  42k—  rj- erf 
B°  So 


( 


2420^ 


(13) 


In  equation  (13),  the  error  function  involving  b  has  been  approximated  for  small  values  of  its  argu¬ 
ment  and  equation  (11)  has  been  used  to  eliminate  the  semi-major  axis  of  the  tolerance  ellipse. 


14 


X  SLOPE  (NORMAL  TO  LINE  OF  SIGHT) 


Figure  5.  A  sequence  of  tolerance  ellipses  for  the  rising  (or  setting)  Sun.  The  receiver  is  looking 
directly  toward  the  Sun  from  an  elevation  of  1  mrad.  The  parameter  next  to  each  ellipse  is  the 
elevation  of  the  solar  center  in  terms  of  the  solar  radius,  e.  The  dashed  line  shows  the  approximation 
for  the  semi-major  axis  of  the  ellipse  given  by  equation  (11). 
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6.  FINAL  RESULT 


Combining  equation  (13)  with  the  result  given  in  appendix  A  for  the  quantity,  BJA,  I  obtain  the 
following  expression  for  the  grazing  reflectivity  of  a  circular  source: 


_T  =  £_  y  1 

r,  ~p  4  a  a/(> f,«) 


2V2  CT  ^ 


(14) 


/(x)=l  +  , 


X2  7  4 

XH - X  +... 

2  24 


Figures  6  and  7  present  a  comparison  of  equation  (14)  with  the  numerical  evaluation  of  equations  (1) 
and  (2)  for  wind  speeds  of  10  and  Va  m  s'1,  respectively. 
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SOURCE  ELEVATION  %  (mrad) 

Figure  6.  Rough  ocean  reflectivity  as  a  function  of  solar  elevation.  The  receiver  is  elevated  by 
1  mrad  and  the  wind  speed  is  10  m  s'1.  The  solid  line  labeled  “Exact”  is  a  numerical  calculation  based 
on  equations  (1)  and  (2).  The  dashed  line  labeled  “Rough”  is  equation  (14).  The  dashed  line  labeled 
“Flat-Top”  is  equation  (23)  of  Zeisse  (1995).  The  dashed  line  labeled  “Cox-Munk/50”  is  equation  (9)  of 
Cox  and  Munk  (1 954)  after  division  by  a  factor  of  50. 
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REFLECTIVITY  RATIO  p 


SOURCE  ELEVATION  Ts  (mrad) 

Figure  7.  The  same  as  figure  6  for  a  wind  speed  of  14  m  s'1. 
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APPENDIX  A 

GRAZING  AREA  OF  A  ROUGH  SURFACE 


The  fractional  area  of  the  rough  ocean  surface  projected  toward  the  receiver  at  any  elevation  is 

:y  =  JJ  (cosCQ/cosen)pd£xdt;y  .  (Al) 

^  co<n  2 

Ur -const 

The  term  in  parentheses  in  the  integrand  of  equation  (Al)  is 

costo/cos  Qn  =Ur-U„/cos  6n 

=  (ara„  +KK  +CrCn)/Cn 

=  ~ar£x  ~brby  +Cr 

=  cos'F^y  +  sinxFr 

where  Ur  and  U„  point  toward  the  receiver  and  along  the  facet  normal,  respectively,  and  ar  vanishes 
because  the  receiver  is  in  the  Y-Z  plane. 

The  behavior  of  equation  (Al)  depends  on  whether  or  not  the  receiver  observes  from  a  grazing 
angle.  In  each  case  the  key  factor  is  the  limit  on  the  angle  of  incidence,  co. 

Well  away  from  grazing  (for  receiver  elevations  above  about  10°),  and  for  typical  wind  speeds,  the 
limit  on  co  may  be  ignored.  That  is  because  under  those  conditions  only  the  front  sides  of  the  facets 
will  be  presented  to  the  receiver  in  the  first  place.  After  inserting  (A2)  into  (Al),  the  trigonometric 
factors  in  (A2)  may  be  brought  outside  the  integral  since  they  do  not  depend  on  facet  slope.  Using  the 
facts  that  (1)  p  is  even,  (2)  the  slopes  are  odd,  and  (3)  p  is  normalized  to  unity,  it  can  be  shown  that 
the  first  term  in  (A2)  vanishes  during  integration  in  (Al)  and  the  second  term  in  (A2)  survives  intact. 
Then  the  integral  in  (Al)  is  equal  to  sin  4V,  the  value  used  by  Cox  and  Munk. 

At  grazing  (for  receiver  elevations  below  about  10°),  the  limit  on  on  may  not  be  ignored. 
Furthermore,  that  limit,  co  =  n/2 ,  is  equivalent  to  the  condition  =  -  tan  xFr .  Introducing  these 

expressions  and  equation  (A2)  into  equation  (Al),  and  using  the  variables  given  in  equation  (10),  the 
integral  in  (Al)  becomes 

1  o°  oo 

—  =  —  Jexp(-M2/2)o/«  J (cosher  v  +  sin 'Fr  )exp(- v' 2 /l)dv,  (A3) 

A  27r  —  —  v„ 

where  the  absolute  value  of  the  reduced  Y-slope  limit  is 

v0  -  tan  vFr  /  o  >  0 .  (A4) 

This  limit  is  small  compared  to  one;  it  has  a  maximum  of  about  1/5  for  the  largest  receiver  angle  of 
10  mrad  and  the  small  wind  speed  of  Va  m  s"1.  After  carrying  out  the  integration  over  u  in  equation 
(A3),  the  projected  rough  area  at  grazing  becomes 


A-1 


"  Ls^7+5£^4 

•nit  {  &  ] 

OO  2 

jvexp(-v2/2)dv  ~  1— — — 
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J=  jexp(-v2/2)<iv  =  JkJi  \  +  erf  (v0 /V2)} 

“V0 


(A5) 


Finally,  I  replace  the  trigonometric  terms  in  equation  (A5)  with  their  small  angle  approximations, 
multiply,  collect  terms,  and  obtain  the  following  expression  for  the  fractional  area  of  the  ocean 
surface  seen  at  the  grazing  angle, 


(A6) 
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